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Analytic solutions are presented for the hydrodynamic collimation of a relativistic fireball by a 
surrounding baryonic wind emanating from a torus. The opening angle is shown to be the ratio of 
the power output of the inner fireball to that of the exterior baryonic wind. The gamma ray burst 
990123 might thus be interpreted as a baryon-pure jet with an energy output of order 10 50 erg or 
less, collimated by a baryonic wind from a torus with an energy output of order 10 52 erg. 
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Recent observations of GRB's and their afterglows appear to indicate that GRB ejecta are collimated. In particular, 
i) the high fluences exhibited by several bursts, notably GRB 990123 for which the fluence indicated an isotropic 
equivalent energy in excess of M Q c 2 |Q, require large focusing factors if the central engine is indeed associated with 
a stellar mass object as widely believed, and ii) the achromatic break seen in the light curve of afterglow emission in 
GRB 990510 Q, can be interpreted (e.g., refs HQ), as due to a transition whereby the blast wave decelerates from 
r > l/ip to r < where tp is the opening angle of the outflow. A blast wave model yielded an opening angle of 
ijj = 0.08 and a corresponding beaming factor of 300 for this source Q. 

If GRB's are indeed collimated, as they appear to be, this would not be surprising. Collimation is a generic feature 
of outflows in astrophysics, and appears in outflows from protostars, AGN, and accretion disks around neutron stars 
(e.g. SS433). There is as yet no universally accepted explanation for such collimation. One possibility is collimation 
by magnetic hoop stress that arises when the field is wound up by rotation. Such a field configuration, however, would 
be kink unstable Q , and in any case is very slowly collimated || . Another possibility is simply that the inertial forces 
in the outer part of the flow collimate material that impacts upon it from within. This hypothesis merely makes 
the plausible assumption that the flow is unsteady or multicomponent, so that there are glancing collisions between 
fluid elements on different streamlines. Even a steady, one component flow can display a collimated jet if it originates 
from a ring, with some of it directed towards the ring's axis of cylindrical symmetry . If that part of the fluid that 
converges at the axis makes a dissipative "sticky" collision, then it remains well collimated as seen by the external 
observer. 

In GRB's, there is an a priori reason to suspect multicomponent flows: If the central engine is a black hole, then 
there needs to be an accretion disk feeding it, and the luminosity from the accretion disk would be so super-Eddington 
that a baryon outflow would be dragged out with the photons and electron-positron pairs [§,U . This would obscure 
the gamma ray burst, so we probably have to say that the gamma ray burst is seen as such only when we look 
down the axis of the baryon poor jet (BPJ) that emerges along field lines connected to the event horizon, which are 
protected by the horizon from baryon contamination ||. The energy that powers the BPJ is either extracted from 
the rotation of the black hole, or deposited there by neutral particles that have crossed over from field lines (such as 
neutrinos that annihilate on the horizon-threading field lines). 

The structure of colliding outflows - We consider an accelerating, baryon poor outflow (BPJ) confined by the pressure 
and inertia of a baryon rich wind. We envision that the baryonic wind is expelled (e.g., from an extended disk or 
torus) over a range of scales much larger than the characteristic size of the central engine ejecting the BPJ, so that 
not too far out its streamlines may diverge more slowly than the BPJ streamlines, thereby giving rise to a collision of 
the two outflows. In general, the collision of the two fluids will lead to the formation of a contact discontinuity across 
which the total pressure is continuous, and two oblique shocks, one in each fluid, across which the streamlines of the 
colliding (unshocked) fluids are deflected. The details of this structure will depend, quite generally, on the parameters 
of the two outflows and on the boundary conditions. 

We seek to determine the structure of the shocked layers assuming that the parameters of the unshocked fluids 
are given. The problem is then characterized by 11 independent variables: the number density, energy density, 
pressure, and velocity of each of the two shocked layers, and the cross-sectional radii of the two shocks and the 
contact discontinuity surface. In what follows subscripts j and b refer to quantities in the BPJ and in the baryonic 
wind, respectively, and subscript s denotes quantities in the shocked layers. To simplify the analysis we approximate 
the shocked layers as cylindrically symmetric, coaxial, one dimensional flows along the channels, and denote by a c , 
a,j, and Of,, respectively, the cross sectional radii of the contact discontinuity surface, the inner shock surface (in the 
BPJ), and the outer shock surface. The energy flux incident into the shocked BPJ layer through the shock is given 
by 

lf k n k = (p j +p j <i , )T j U j Bm5 j , (1) 

and the momentum flux in the axial direction incident through the shock is 

Tfm = pj + (pj + pj c 2 ) Uf cos 8j sin 5j. (2) 

Here rik are the components of a unit vector normal to the shock surface, 0j is the angle between the velocity of the 
fluid just upstream the shock and the BPJ axis, Tj and Uj are the Lorentz factor and 4-velocity of the unshocked 
BPJ, and Sj is the deflection angle (i.e., the angle between the velocity of upstream fluid and the shock surface), 
which is related to 6j through: sinSj = (sin 6^ + cos 9 jdaj/dz)/^/i + (daj/dz) 2 . With the above results the energy 
and momentum equations of the shocked fluid can be written as 

^-[T^nia 2 - a 2 )} = T° k n k dS/dz = ( Pj + Pj c 2 )T 'jUj-n a j (sin 6 'j + cosOjddj/dz), (3) 
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and 



^-[T 7 --"K{a 2 c - a 2 )] = TfmdS/dz = pj + (pj + PjC 2 )U 2 cos 0,- to.,- (sin 0,- + cos 'j da j/dz), (4) 

where X^j 2 , are the energy and momentum fluxes in the axial (z) direction of the shocked BP J fluid, and dS = 
najWl + (daj /dz) 2 dz is an element of the shock surface between z and z + dz. The projection of the momentum 
equation on the direction normal to the shock surface yields 

p sj = rijUk = Pj + {pj + PjC 2 )Uf sin 2 5j. (5) 

The last equation holds provided the transverse momentum flux incident through the shock exceeds the free expansion 
pressure of the shocked fluid. The above equations must be supplemented by an equation of state, p S j — Psj(psj), 
and a continuity equation that governs the density profile of the shocked fluid. Likewise, the equations describing the 
structure of the shocked baryonic wind can be written in the same form with the subscript j replaced by b. Finally, 
the requirement that the net momentum flux across the contact discontinuity surface vanishes implies that p s b — Psj ■ 
A complete treatment requires numerical integration of the coupled set of 11 equations derived above, subject to 
appropriate boundary conditions, and is left for future investigation. Below, we consider for illustration a simple 
situation where the BPJ is assumed to be pressure dominated, and the assumption is then justified a posteriori. This 
corresponds to the limit aj = smSj = 0. We further suppose that the shocked baryonic layer is very thin, so that to a 
good approximation we can set a b = a c in the above equations. Adopting a relativistic equation of state for the BPJ, 
p = pc 2 /3, mass, energy and momentum conservation imply: 

p:/y>-«; = c, (6) 

(p s j + p S jC 2 )r 2 sj ira 2 c c ~ 4p sj T 2 sj wa 2 c c = Lj, (7) 
Psj =Pb + (Pb + PbC 2 )U 2 sin 2 S b , (8) 

where C and Lj, the corresponding BPJ power, are constants. Equations (|6|), (Q) can be combined to yield a relation 
between the Lorentz factor, the pressure, and the cross sectional radius: 

(T s j/T sjo ) = (a c /a co ) = [psj/PsjoY 1 ^, (9) 

where T S j a , a co , and p S jo are the corresponding values at z = z Q . On substituting eq. @ into eq. (||), one obtains a 
differential equation for a c (z): 

4 6 + Pb<?)Ul (sin 6 + cos9 b da c /dz) 2 

{ac/aco) =Pb Psjo + r - , 21 • (10) 

Psjo [l + {da c /dz) 2 \ 

Pressure confinement - We consider first the possibility that the BPJ is collimated by the pressure of the baryonic 
wind. Such collimation will occur naturally in the region located within the acceleration zone of the BPJ and above 
the acceleration zone of the baryonic wind (roughly above its critical point) , even if the two outflows have the same 
equation of state. The reason is that in this region the density profile of the BPJ declines with radius as r~ 3 (for a 
conical BPJ) while that of the wind declines as r~ 2 . Taking 5b = in equation m) and using eq. (||) one finds, 

a c {z) = a co [p b (z)/p bo }~ 1/4 . (11) 

Adopting an equation of state of the form p b {n b ) oc for the baryonic fluid yields, a c {z) — a co [n b {z)/n ob ]~ a / i . In 
case of a conical (or spherical) wind, n b (z) oc z -2 , and the cross sectional radius scales as a c {z) oc z Q / 2 . Convergence 
occurs for a < 2, which is the case for both relativistic gas (a = 4/3) and non-relativistic gas {a = 5/3). 

The above analysis does not take into account effects associated with the formation of a rarefaction wave in the 
baryonic outflow due to the divergence of its stream lines near the interface separating the two fluids. This will 
result in a steeper decline of the pressure supporting the BPJ at radii where the wind is highly supersonic, and the 
consequent alteration of the profile of the contact discontinuity surface. We naively expect collimation to take place 
predominantly over a range of radii at which the Mach number of the baryonic wind is mild, unless some fraction 
of the wind power is dissipated, e.g., due to formation of shocks, far enough out. This suggests that only a modest 
beaming factor may be attainable in this case. The determination of the level of collimation in this scenario is left for 
future work. 
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BP J confinement by a wind emanating from a torus - As a second example we consider a baryonic wind emanating 
from a thin torus of radius R located at z = 0, where z is the axial coordinate, and centered around the central engine 
ejecting the BP J (see Fig. 1). We suppose that on scales of interest the baryonic wind is highly supersonic, so that 
momentum transfer from the wind into the BPJ is dominated by ram pressure. We can therefore neglect the kinetic 
pressure, p b , m cq. ©. Now the ram pressure of the baryonic outflow just upstream the oblique shock at some 
position z is related to the total wind power , L b , through: 

(p b + Pb<?)Ul = (3 b L b /{4ir 2 ca c r), (12) 

where r = [z 2 + (R — a c ) 2 ] 1//2 is the distance between a point on the torus and the nearest point to it at the shock 
(see figure 1), and (3 b is the (terminal) wind velocity. The angle between the wind velocity and the BPJ axis is given 
by tan#b = (R. — a c )/z. Substituting eq. ( p"2"|) into eq. (|To|) , and using eq. (fjj) yields, 

{R - a c + da c /d In z) 2 _ 2 _ t [z 2 + (R - a c ) 2 ] 3 / 2 

l + (da c /dz) 2 ° coX al [6) 

with x — 7r ^ 1 Pb^ 2 j {Lb/ Lj) being approximately the ratio of baryonic wind and BPJ powers. At large enough axial 
distances from the plane of the torus, z >> R, an approximate solution to eq. (^) can be obtained by expanding 
the various terms in powers of R/z, and becomes more exact at large z. To second order the BPJ in this regime is 
conical, viz., a c (z) = b + ctz, where a, the opening angle, and b are constants, and satisfy the relation: 

X (R/a co - b/a co ) 2 = 1 + 3 > . (14) 

or 

The parameter b/a co depends on the properties of the solution near the central engine (at z < R) and must be 
determined numerically. The requirement that the BPJ is confined by the ram pressure of the baryonic outflow at 
z = yields, 

(R/a co -1)= X , (15) 
where it has been assumed that da c (z = 0)/dz — 0. Combining the last two equations, one finds 

1 J = X {x + l-b/a co ) 2 . (16) 
For x >> 1 the latter equation implies that 

a-x^ + Oix- 2 )- (17) 

Numerical integration of eq. ( |l3| ) confirmed that the BPJ is indeed conical with an opening angle given by eq. (|l^) , 
except for a small range (z < R) near the central engine. BPJ profiles computed numerically for different values of 
the parameter x are exhibited in Fig. 2, and the corresponding dependence of a on x is displayed in Fig. 3. 

Once the BPJ Lorentz factor exceeds a" 1 ~ x> it will remain conical with the same opening angle regardless of 
the external conditions. This occurs at a distance z = i?x 2 /(l + x) from the central engine. Thus the baryonic wind 
must extend to at least this scale. 

The apparent luminosity, L app , exhibited by a beamed source having a power Lj is roughly L app ~ f~ l Lj, where 
/ = AS1/2-7T is the corresponding beaming factor. For the above model we find 

/ = 7r a 2 /27r = 1/2 X 2 - (ir 2 /2)(L 3 /L b ) 2 , (18) 

and 

L app ~ {2/n 2 )(L b /L ] )L b . (19) 

In conclusion, we have shown [equation (|l9|)] that an inner relativistic jet collimated by an outer modestly relativistic 
jet can appear brighter to the observer within the beam by a factor that is inversely proportional to the intrinsic 
power of the inner jet. When all else is the same, a low luminosity inner jet actually appears brighter than a high 
luminosity one to the observer within its beam. The cost, of course, is that, fewer observers are within the beam. 
Ironically, most of the energy of the GRB is the unseen outer jet. 
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Now consider the remarkable GRB 990123, which had an isotropic equivalent energy output of about 3 x 10 ergs. 
Applying equations through ([l9]), we propose that the actual energy in the gamma ray fireball was only about 
77 2 10 50 ergs, (it may have even been smaller than average) and an outer baryonic wind output of order 77IO 52 5 erg. 
The corresponding beaming factor is / = 10~ D 7r 2 77 2 /2. Here r\ is assumed to be of order unity or less, otherwise the 
energetics strain theoretical models. 

A rough observational lower limit on 77 can in principle be derived from the consideration that GRB 990123 could 
not have a collimation angle that is implausibly smaller than for average GRB's. Although it was an unusual burst, 
the a priori chance of seeing it is unlikely to have been less than 10 -3 , and, considering that it is a BcppoSax-located 
GRB and that there has been at least one other burst of comparable (if somewhat smaller) intrinsic magnitude, 
probably considerably more than 10 -3 . The largest uncertainty is the beaming factor for the average GRB. If we take 
that as being of order 10~ 3 , corresponding to an opening angle of about 0.06 radians, this suggests a value for 77 of 
at least of order 0.3. This means that the minimum energy requirement for GRB 990123 is still of order 10 52 ergs 
or more. On the other hand, the fact that it can come out as a baryonic wind relaxes the demands on theoretical 
models. For example, a magnetocentrifugal wind could plausibly release several percent of a solar rest energy over 30 
seconds, whereas it would be hard to get this much from v — v annihilation. 

Another example is GRB990510, for which the achromatic steepening of the optical and radio light curves observed 
after about 1 day ||, can be modeled as due to the evolution of a jet with an opening angle of ~ 0.08 rad, and 
a corresponding beaming factor of 300. The prompt GRB emission should be at least as beamed as the afterglow 
emission, but could in fact be more beamed. The isotropic gamma-ray emission inferred for this source is 3 x 10 53 
ergs. Employing equations (|l8|) and ([l9|) yields a total energy of < 10 51 ergs for the BPJ and 3 x 10 52 ergs for the 
confining wind. 

We acknowledge support from the Israel Science Foundation 



Galama, T.J., et al. Nature, 398, 394 (1999) 
Harrison, F, et al. Astroph. J., 523, L121 (1999) 
Rhoads, J.E. Astroph. J, 487, LI (1997) 



Moderski, R., Sikora, M., & Bulik, T. 1999, ]astro-ph/990431(\ (1999) 

Eichler, D. Astroph. J, 419, 111 (1993) 

Begelman et al. Astroph. J, 505, 835 (1998) 

Peter, W. & Eichler, D. Astroph. J, 466, 840 (1996) 

Levinson, A., & Eichler, D. Astroph. J, 418, 386 (1993) 

Woosley, S.E. Astroph. J, 405, 273 (1993) 



4 



FIG 1 




FIG. 1. Schematic illustration of BP J confinement by the ram pressure of a baryon rich wind (represented by the dashed 
lines) emanating from a torus of radius R. The torus is centered around the central engine ejecting the BP J. Fluid elements 
of the baryonic outflow collide with the BPJ at an angle S(z), and are deflected when passing through an oblique shock. The 
shocked wind layer is assumed to be very thin. At distances larger than roughly R the BPJ becomes conical with a semi-opening 
angle a. 
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FIG. 2 
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FIG. 2. Cross sectional radius of the BPJ versus axial distance from the central engine, computed numerically for different 
values of X- 
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FIG. 3 
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FIG. 3. A plot of q 1 against x, obtained from the numerical integration of eq. 
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